Search for: All records

Marine heatwaves (MHWs) are warm sea surface temperature (SST) anomalies with substantial ecological and economic consequences. Observations of MHWs are based on relatively short instrumental records, which limit the ability to forecast these events on decadal and longer timescales. Paleoclimate reconstructions can extend the observational record and help to evaluate model performance under near future conditions, but paleo-MHW reconstructions have received little attention, primarily because marine sediments lack the temporal resolution to record short-lived events. Individual foraminifera analysis (IFA) of paleotemperature proxies presents an intriguing opportunity to reconstruct past MHW variability if strong relationships exist between SST distributions and MHW metrics. Here, we describe a method to test this idea by systematically evaluating relationships between MHW metrics and SST distributions that mimic IFA data using a 2000-member linear inverse model (LIM) ensemble. Our approach is adaptable and allows users to define MHWs based on multiple duration and intensity thresholds and to model seasonal biases in five different foraminifera species. It also allows uncertainty in MHW reconstructions to be calculated for a given number of IFA measurements. An example application of our method at 12 north Pacific locations suggests that the cumulative intensity of short-duration, low-intensity MHWs is the strongest target for reconstruction, but that the error on reconstructions will rely heavily on sedimentation rate and the number of foraminifera analyzed. This is evident when a robust transfer function is applied to new core-top oxygen isotope data from 37 individualGlobigerina bulloidesat a site with typical marine sedimentation rates. In this example application, paleo-MHW reconstructions have large uncertainties that hamper comparisons to observational data. However, additional tests demonstrate that our approach has considerable potential to reconstruct past MHW variability at high sedimentation rate sites where hundreds of foraminifera can be analyzed. 

Abstract Scientific understanding of low-frequency tropical Pacific variability, especially responses to perturbations in radiative forcing, suffers from short observational records, sparse proxy networks, and bias in model simulations. Here, we combine the strengths of proxies and models through coral-based paleoclimate data assimilation. We combine coral archives ( δ 18 O, Sr/Ca) with the dynamics, spatial teleconnections, and intervariable relationships of the CMIP5/PMIP3 Past1000 experiments using the Last Millennium Reanalysis data assimilation framework. This analysis creates skillful reconstructions of tropical Pacific temperatures over the observational era. However, during the period of intense volcanism in the early nineteenth century, southwestern Pacific corals produce El Niño–Southern Oscillation (ENSO) reconstructions that are of opposite sign from those from eastern Pacific corals and tree ring records. We systematically evaluate the source of this discrepancy using 1) single-proxy experiments, 2) varied proxy system models (PSMs), and 3) diverse covariance patterns from the Past1000 simulations. We find that individual proxy records and coral PSMs do not significantly contribute to the discrepancy. However, following major eruptions, the southwestern Pacific corals locally record more persistent cold anomalies than found in the Past1000 experiments and canonical ENSO teleconnections to the southwest Pacific strongly control the reconstruction response. Furthermore, using covariance patterns independent of ENSO yields reconstructions consistent with coral archives across the Pacific. These results show that model bias can strongly affect how proxy information is processed in paleoclimate data assimilation. As we illustrate here, model bias influences the magnitude and persistence of the response of the tropical Pacific to volcanic eruptions.